AD-A050  SaS  CORNELL  UNJV  XTHACA  N Y SCHOOL  Op  ELECTRICAL  EnOINEERINO  F/g  18/1 

preliminary  study  of  HISh  neutron  flux  fusion  HEaTIN8.(U) 

*<icfctrn  77  R L LIBOFF  AFOSR-77-3326 

UNCLASSIFIED  AFOSR-TR-78>0257  NL 


AD 

A060886 


4-78 


4.  monitoring  agency  name  • ADDRESr^// d///eren»  Irom  Controlling  Ollice)  I 15.  SECURITY  CfASS.'fo/fcifS-report^ 


AFOSR/NP 

Building  UlO,  Bolling  AFB 
Washington,  D.C.  20332 


16.  DISTRIBUTION  STATEMENT  (ol  Ihit  Rtporl) 


Approved  for  public  release; 
distribution  ualimited. 


UNCLASaiFIEy  n 


ISa.  DLCLASSIFICATION/ downgrading 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (oi  the  ebatrect  entered  In  Block  20,  U different  from, Report) 


19.  KEY  WORDS  (Continue  on  reverse  aide  it  neceaaary  and  identify  by  block  number) 


Neutrons,  Fusion,  Heating,  Lattice-deformation,  Collision-loss. 


i^O.  abstract  (Continue  on  reverse  aide  If  neceaaary  and  Identify  by  block  number)  ,, 

A heating  scheme  for  nuclear  fusion  is  proposed  based  on  the  availability  of  a 
high  flux,  low  energy  neutron  source.  The  heat  is  derived  in  the  reaction 
L#\n^)He'^  resulting  from  the  incidence  of  a low  energy  neutron  beam  on  a sampld 
of  Li^^.  The  energy  release  per  reaction,  Q = U.6  MeV,  is  converted  through  a 
electron  Coulomb  collisions  thereby  quickly  dissociating  the  solid  sample  to  the 
plasma  state.  For  ^ eV  neutrons  it  is  estimated  that  this  dissociation 

occurs  in(,r  7 msec  for  an  incident  flux  of  ^ 10^  cm~^  sec~^^.  The  possibility  oi 
furtlier  dpi  v.i ng  the  heal, oil  fuel  t.n  fnsiop  is  also  di  reussed . ^ 


DD  1473  EDITION  OF  I NOV  65  IS  OBSOL^E 


t bgo'f 


CthttJ’  'Ibfo'U  p 


C »•  / 5^ 

llNrT.ASr^lFTFD 


November  3,  1977 

Report  No.  1-77 
AFOSR 


Approved  for  public  release; 
dlstrlbutloa  unlimited* 


4 


AIR  rORCX  OFFICE  OF  SCIENTIFIC  RESEARCH  (AFSC) 
NOTICE  OF  TRANSMITTAL  TO  DDC 
This  technical  report  has  b-oa  rcvlcwoJ  end  Is 
approved  Tor  pitbl:C  rclcjiso  IA'»/  Ai^R 
Distribution  is  unlimited. 

A.  0.  BLOSE  , . 

Taohnloal  Infonatlon  Officer u-b-cu;  . - • 


/ 


Abstract 


A heating  scheme  for  nuclear  fusion  is  proposed  based  on  the 

availability  of  a high  flux,  low  energy  neutron  source.  The 

heat  is  derived  in  the  reaction  Li  (n,T)He  resulting  from  the 

incidence  of  a low  energy  neutron  beam  on  a sample  of  Li^D. 

The  energy  release  per  reaction,  Q = 4.6  MeV,  is  converted 

through  electron  Coulomb  collisions  thereby  quickly  dissociating 

-3 

the  solid  sample  to  the  plasma  state.  For = 10  eV  neutrons 

it  is  estimated  that  this  dissociation  occurs  in  = 7 msec  for  an 
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incident  flux  of  = 10  cm  sec  . The  possibility  of  further 
driving  the  heated  fuel  to  fusion  is  also  discussed. 
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1.  Introduction 

Recent  world  wide  efforts  have  been  initiated  to  obtain  neutron 

• IT  2 

fluxes  of  the  order  of  10  neutrons/cm  -sec,  employing  high  energy 

pulsed  proton  bombardment  of  heavy  atoms.  Present  day  thermal 
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neutron  fluxes  of  about  10  neutrons/cm  -sec  are  generated  by 

nuclear  reactors  at  Brookhaven,  Grenoble  and  Oak  Ridge  [1,2]. 

In  this  paper  we  propose  a scheme  for  fusion  heating  which  is 

based  on  the  availability  of  high  neutron  fluxes  of  low  energy.  In 

the  envisaged  scheme  a beam  of  low  energy  neutrons  is  incident  on  a 

sample  of  Li^D.  The  low  energy  neutrons  are  readily  captured  by 
0 

the  Li  nuclei  and  release  an  a particle  and  a triton  through  the 
reaction 

(1)  Li^(n,T)He‘* 

The  Q value  of  this  reaction  is  4.6  MeV.  The  T particle  emerges  with 
2.63  MeV  and  the  a particle  with  1.97  MeV.  This  energy  goes  into 
heating  the  sample  and  ionization.  In  bringing  the  sample  to 
fusion  temperature,  it  is  necessary  first  to  exhaust  the  crystal 
energy,  thereby  effecting  a fusion  plasma  comprised  of  ionized  Li, 

D,  He  and  T elements. 

2.  Collision  Loss 

The  emitted  particles  in  reaction  (1)  lose  energy  primarily 
to  electron  collisions  and  ionization  [3].  The  collision  loss  per 
unit  length  is  given  by  the  Bethe-Bloch  formula  [4] 


^wnppwimRppii 


T 


(2) 
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dx  " eE 


z*  In  gE 


g = 4e/ZI,  e = m/M 

I = 11.5  eV 

where  z is  the  charge  number  of  the  incident  particle.  The 
ionization  loss  is  given  by 


(3) 
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dx  2 In  gE  dx 


With  electron  density  n^  = 2.37  x 10^^  cm"^  and  atomic  number  Z 
we  obtain  for  the  triton  loss. 
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dx 


e _ 170 
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(MeV/cm) 


Let  us  assume  a mean  triton  energy  of  1.32  MeV  relevant  to  its 
slowing  down  trajectory.  Then  (4)  gives,  for  the  triton 


(5) 


dx 


dEj 

dx 


428  MeV/cm 


= 0.15 


dx 


With  an  assumed  mean  energy  of  0.985  MeV  for  the  a particle  we 
obtain 
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d£ 


dx 


C _ 907 


In  16  E = 2540  HeV/cm 


(6) 


dE_  dEp 

= 0.18 

dx  dx  . , . _ ^ 

Radiation  loss  may  be  neglected  compared  to  collisional 

loss  [5]  providing  E/MC  <<  1.  This  inequality  is  obeyed 

for  the  case  at  hand.  This  observation  together  with  (5) 

and  (6)  indicate  that  we  may  assiime  the  primary  loss  mechanism 

to  be  heating  of  the  lattice  through  ion-electron  collisions. 

3.  Lattice  Deformation  Time 


The  lattice  energy  of  LiD  is  [6 ] 

(7)  Ej^  = 218.8  Kcal/mole 

= 5.71  X 10^®  MeV/mole 

The  completion  of  this  energy  leaves  the  fuel  sample  in  a state  of 
isolated  Li^  and  D“  ions.  With  this  energy  consumed  collisional 
decrement  is  still  the  primary  means  of  energy  loss , which  now  may 
be  assumed  to  go  into  heating  the  sample . 

If  the  original  sample  thickness  is  much  less  than  the  neutron 
mean  free  path  in  reaction  (1),  we  may  assume  that  the  incident 
neutron  flux  is  roughly  maintained  through  the  sample.  In  this  case 
the  neutron  reaction  rate  is 

(8)  R = JA£na  = JnaV(sec”^) 

I 

\ 

' where  n is  Li  number  density  and  V is  the  volume  of  the  LiD  crystal 

! which  has  cross  secional  area  A and  thickness  Z.  The  cross  section 


4,6  Mey,  which  with  (8)  gives  the  yield 

3 

y = 4.6  Jno  MeV/cm  -sec 

3 

The  molar  volume  of  LiD  is  10.1  cm  /mole  [6]  and  the  last  equation 
gives  the  molar  yield 

(9)  Y = 46  Jno  MeV/mole-sec 

Comparison  with  (8)  indicates  that  the  lattice  energy  is  supplied 
in  the  time 

(10)  t,  = ^ 

Under  the  somewhat  severe  assumption  that  a means  of  confinement 
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can  be  applied  to  roughly  maintain  the  fuel  size,  then  n = 6 x 10 

_3 

cm  and  the  last  equation  becomes 


2 X 10' 


(sec) 


From  low  energy  neutron  data  [7 ] relevant  to  reaction  (1)  one  may 
construct  the  expression 

C12)  o = 135  X (E  in  eV,  o in  barns) 

_3 

For  incident  neutron  energy  ~10  eV  this  latter  formula  gives 

3 

o = 3 X 10  b which  with  (11)  gives  the  lattice  deformation  time 
(13)  t^  = 

It  follows  that  for  J = 10^^,  t^  = 7 msec.  Once  the  lattice  heats, 
the  relative  n-Li  kinetic  energy  increases  thereby  decreasing  a 
and  increasing  t^^ . 

4.  Fusion  Time 

In  the  completely  ionized  state,  each  LiD  molecule  contributes 


m 
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six  particles  to  the  plasma.  Thus  one  mole  of  original  fuel  sample 

24 

is  converted  to  3.6  x 10  particles  in  the  completely  ionized 
state.  The  T-d  fusion  reaction  has  its  peak  cross  section  [8]  at 
0.1  MeV.  At  thermal  equilibrium  this  corresponds  to  the  fusion 
molar  energy 


= 3.6  X 10^^  MeV/mole  = 6.3  x 10** 


With  (9),  we  then  obtain  the  fusion  time 

E. 


(14) 


'f  _ 0.13 
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At  E = 10  eV,  this  gives  the  time 

.19 


(15) 


_ 4.3  X 10 


so  that  t^  or  1 sec  for  J o'  10 


X 
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( sec) 


This  magnitude  is  somewhat  beyond 
the  borderline  of  present  day  neutron  flux  values  . Once  again  we  note 
that  this  time  estimate  is  low  due  to  the  increase  in  the  relative  n-Li 
kinetic  energy  at  elevated  temperature. 

5.  Sample  Thickness  and  Conclusion 


The  neutron  mean  free  path  for  the  reaction  (1)  is 

,-23 


(16) 


0 = 1_  = i-7  X 10 
*'n~  no  ~ a 
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At  E = 10  eV  (corresponding  to  velocity  of  400  m/sec)  the 
reaction  mean  free  path  is 


£ - 60  pm 
n 


As  stated  previously  this  is  the  approximate  width  of  the  original 
fuel  sample,  which  must  be  large  compared  to  the  thermalization 
length 
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From  (5)  and  (6)  we  see  that  the  triton  has  the  longer  thermalization 
length.  At  triton  energy  of  1.3  MeV  we  obtain 


= 30  pm 

which  is  approximately  half  the  width  of  the  sample.  At  larger 
neutron  energy,  this  inequality  becomes  more  favorable,  however  both 
tj^  (11)  and  tj  (14)  also  increase  with  larger  neutron  energy 
rendering  the  process  still  more  unfeasible. 

In  summary,  our  rough  estimates  have  led  to  the  following 

_3 

conclusion.  Namely,  that  a beam  of  10  eV  cryogenic  neutrons 

incident  on  a thin  sample  of  Li^D  of  thickness  2i60  ym  will  act  as 

a catalyst  for  nuclear  reactions  which  will  in  turn  heat  the  sample. 

At  these  values,  thermonuclear  temperatures  may  be  reached  in  ^ 1 

19  -2 

second,  provided  a neutron  pulse  with  flux  in  excess  of  10  cm 
sec  ^ and  means  of  confinement  are  available.  In  the  absence  of  this 
confining  mechanism,  it  may  still  prove  feasible  to  employ  the 
initial  heating  scheme  in  conjunction  with  more  pedestrain  plasma 
confinement-compression  devices  [9,10]. 

A mesure  of  the  feasibility  of  this  approach  is  offered  by  the 
Lawson  criterion  [ll].  This  criterion  stipulates  that  the  product 
nx  of  particle  density  n and  confinement  time  x for  the  D-T  reaction 
must  be  in  excess  of  lO'*'  s-cm~  for  positive  energy  gain.  Present 
day  fusion  devices  have  not  attained  this  value.  In  the  heating 

process  of  the  present  scheme  we  ask  that  the  fuel  element,  at 

23  -3  20 

n ~ 10  , be  maintained  for  ~ 10  sec.  which  gives  nx  = 10  . This 

large  magnitude  expresses  the  difficulties  inherent  to  the  proposed 
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